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We compare the approach to the chiral transition line K c (/3) in quenched and full 
compact lattice QED with Wilson fermions within the confinement phase, especially 
in the pseudoscalar sector of the theory. We show that in the strong coupling limit 
(j3 = 0) the quenched theory is a good approximation to the full one, in contrast 
to the case of /3 — 0.8. At the larger /3-value the transition in the full theory is 
inconsistent with the zero-mass limit of the pseudoscalar particle, thus prohibiting 
the definition of a chiral limit. 



x. 

1 Introduction and model description 

Chiral symmetry as a major concept in continuum quantum field theory has re- 
mained a problematic topic in lattice gauge theories over the years. It is well-known, 
that for Wilson fermions chiral symmetry is explicitly broken in QCD and QED on 
the lattice [0, 0. Hopefully, it can be recovered by fine-tuning of parameters in 
the continuum limit. Then some line k c (0) in the phase diagram is associated 
with the chiral limit of the theory. On the other hand, for non-vanishing lattice 
spacing only a partial restoration of chiral symmetry at k — k c (/3) is possible with 
Wilson fermions 0, ||. It is still an open question, how this mechanism of partial 
symmetry restoration should eventually be integrated into the general conception 
of spontaneously broken chiral symmetry. One cannot exclude that the breakdown 
of some other symmetry group governs the dynamics of the transitions at k c (/3) 
(e.g., p|). Viewed in this light the vanishing of the pseudoscalar 'pion' mass m n 
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for k — > k c (/3) is a necessary but not sufficient condition for probing the chiral 
limit. Another point which sharpened the look on the chiral limit in QCD || [jj 
is the discussion of 'enhanced logs' due to quenching, demonstrating the role of 
dynamical fermions in chiral properties of the theory. 

In this letter we are concerned with the behavior of fermionic observables close to 
k c (/3) in the confinement phase of compact QED. We confront the full theory with 
its valence approximation. Though the non-perturbative regime of compact QED 
is itself an interesting subject the close analogy with non-Abelian gauge theories 
makes it also a valuable test ground for QCD with Wilson fermions. 

In the rest of the section we will introduce the model and give the main notations. 
Starting point is the partition function of 4d compact QED which reads as follows 

Zqed = J [dt/][##] e~ s ^ u > , (1) 
where Sw(U, ip, ip) denotes the standard Wilson lattice action 

S W = S G (U) + S F (U,^^) (2) 
consisting of the plaquette action 

S G (U)=(3- J2 (1-cos^) , (3) 

x;fi>u 

and the fermionic part SpiU^^ip) 



f 



f=l x,y 



M. xy = 1 — K 



5y )X+il ■ (1 - 7 M ) • U xli + 6 yiX -fi ■ (1 + 7 M ) • Ul« , (4) 



with (3 = \/gl are) and U x ^ = exp(i9 Xfl ), 6 Xfl £ (— n,ir] represent the link fields. The 
plaquette angles 9 x;fiV in eq.(|) are given by 9 X .^ U = 9^^ + 9 x+fi;u - 9 x+i/ .^ - 9 X . V . 
In the fermionic part of the action Ai xy denotes Wilson's fermionic matrix with 
the hopping parameter k and the flavor-index /. 

The phase diagram of this model has been studied in ||. Within the region of 
values < k < 0.30 the existence of four phases has been established. The line 
k c (/5) separates both the confinement phase ( < (3 < [3q with /? ~ 1.01 ) and 
the Coulomb phase ( (3 > /3 ) from two upper phases (which we called the 3 rd phase 
at weak and the 4 th phase at strong coupling). At /3 = 0.8 it was observed that the 
scalar condensate (4>ip) has a discontinuity at k c , therefore we called the transition 
at this point a 1 st order transition. However, the ^-dependence of the pseudoscalar 
observables was not studied in detail. 

The fermionic observables to be discussed are the pion norm 
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(n> = -L ■ (tt(m- 1 15 m- 1 15 )) g , (5) 

being a good indicator for small eigenvalues of the fermionic matrix and the mass 
of the pseudoscalar particle m n , which is extracted from the pseudoscalar zero- 
momentum correlator 




r(r) = (^h^(r,£) ■ ^75^(0, y)) 

s x,y 

In eqs.([|||) ( )g indicates averaging over gauge field configurations, and V = N T ■ 
Ng is the number of sites. Sp means the trace with respect to the Dirac indices. 
Other observables like (p m0 n) - the density of DeGrand-Toussaint monopoles || 
- and the photon correlator have also been determined. 

Formally we define the bare fermion mass parameter m q by 

For the determination of the below. 

In this work we discuss data from an 8 3 x 16 lattice at two values of (3 within 
the confinement phase. For the production of dynamical gauge field configurations 
we employed a CRAY-T3D implementation of the Hybrid Monte Carlo (HMC) 
method. A detailed presentation of algorithmic issues, like the tuning of the HMC 



parameters when approaching k c is deferred to JlC] 



2 Effects of dynamical Wilson fermions 

First, we shall discuss the behavior of the bulk observable (II) when approaching 
the line k c (/3) at fixed (3 within the range < (3 < j3 (confinement phase). As 
in our previous work for representative /3-values we have chosen (3 = 0.8 and the 
strong coupling limit (3 = 0. In the latter limit the comparison with analytical 
results is possible (e.g., || PHI ). 

Provided the pseudoscalar mass vanishes for m q — > it will yield the dominant 
contribution to the pion norm (II) ~ l/ m ^- I n case °f a PCAC-like relation 
between m n and m„ the pion norm can be expressed in the following form 



(H) = 9l + Cl , m q ^0, (8) 
m q 
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where Co > - up to a factor - is the subtracted chiral condensate [jn| (see also 
fi"3H). C\ = (Il m ) is the contribution of the massive modes. 

In Figure [I] data for the quenched and dynamical cases at (3 = are presented 
to view. The quenched approximation follows the m^-dependence of the full pion 
norm (IT) very well, even quantitatively. As m q — > + the singularity of (II) at 
(3 = is well described in both cases by eq.(H) with values of Co, C\ listed in Table 
|I|. Note that Co and G\ differ somewhat for the quenched and dynamical cases 
and could not be forced to coincide by shifting k c . The quenched and dynamical 
theories nevertheless exhibit the same functional dependence on m q when m q — > 0. 

For k^j k c the averages of fermionic bulk observables in the quenched approx- 
imation become poorly defined due to large fluctuations caused by 'exceptional' 



configurations |13| , 14]. In the theory with dynamical fermions we also observe in- 
creasing fluctuations of, e.g., II when k is tuned towards k c from below. We 
were able to proceed to k = 0.242 {m q ~ 0.0253) at (3 = 0, before being faced 
with serious problems with the acceptance rate of the HMC method. 

The situation changes drastically when the gauge coupling (3 is increased. By 
examining time histories of II and other observables at (3 = 0.8 we observe the 
formation of metastable states in a 'critical' region around k c in the presence of 
dynamical fermions. Figure 0a illustrates a clearly double peaked distribution of 
II close to k c for the case of dynamical fermions. 

As an example concerning the behavior of other observables, the dependence of 
(Pmon) on n is depicted in Figure §b. The evolution of (p m on) resembles the 
situation of the confinement-deconfinement transition at (3q for a sufficiently small 
fixed k and varying (3 [|J. 

Measurements of the effective photon energy extracted from plaquette-plaquette 
correlators for non-zero momentum confirm that in the case of dynamical fermions 
the system undergoes a confinement-deconfinement transition at n c (see 0). With 
increasing k the effective energy of the photon rapidly decreases around k c and 
becomes well consistent with the lattice dispersion relation for a zero-mass photon. 

We used further 'order parameters' to determine k c and to make sure that there 
are no other transition points different from that within the investigated K-range. 
For example, the variance <t 2 (II) which is a suitable parameter to locate the line 
K c (f3) within the Coulomb phase || peaks at the same k c . 

In Figure ^| as counterpart of Figure [I] we confront the dependence of (II) on 
m q in the quenched approximation with the corresponding data when dynamical 
fermions are taken into account. The general features of the valence approximation 
at (3 = 0.8 are the same as at (3 = 0. (II) has a singularity for m q — > + 
described very well by eq.(||) (dashed line in Figure |J) with parameters C ; C\ 
given in Table |T|. However, the effect of the fermionic determinant is now seen as 
a qualitative change in the behavior of the pion norm, which does not behave as 
~ l/m q but rather exhibits a finite discontinuity accompanied by a metastable 
behavior (cf. Figure 0a) around k c . In the quenched theory averages of II become 
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statistically not well-defined for the case of (3 — 0, while in the 

dynamical case at /3 — 0.8 there is no problem to go beyond k c (i.e., into the 



3 phase [ 
around K r 



The dependence of (IT) on k (respectively m q ) is not symmetric 



To substantiate the emerging picture at (3 = and (3 
the evolution of the pseudoscalar mass m n when k —>■ k c . 



= 0.8 we will discuss 
We present in Figure 



|] the dependence of on k for the full and quenched theories on an 8 x 16 
lattice at (3 = 0. The quenched data for k very close to obtained by an 

improved estimator of []T3| in order to increase the signal-to-noise ratio. Since 
the existence of a massless pseudoscalar particle is predicted by strong coupling 



arguments H] we extrapolate our dynamical data for 



linearly to zero in order 
In Table 



to determine k c ((3 = 0) as we have done in the quenched case |L3] . 
[l] we list the values of k c obtained for the dynamical and the quenched cases 
on an 8 3 x 16 lattice. The extrapolated values of k c for the quenched case 
are well consistent with the prediction at strong coupling [[|. In both, quenched 
and dynamical cases we observe the following dependence of m\ on the hopping 
parameter when approaching k c from below 



ml ~ 



K 

1 



which in this limit transforms into a PCAC-like relation between m\ 
fermion mass m n : 



(9) 

and the bare 



B ■ m 



(I ! 



+ 



(10) 

The corresponding slopes B for the quenched and full theories coincide within the 
errorbars (see Table 0). Thus, Figure § suggests a zero- mass pseudoscalar particle to 
exist. However, this is a necessary but not a sufficient prerequisite for the definition 
of the chiral limit. 



p = 






B 


Co 


Ci 


dynamical 


0.2450(6) 


4.87(5) 


0.895(1) 


0.88(1) 


quenched 


0.2502(1) 


4.91(4) 


0.996(2) 


0.80(1) 


(3 = 0.8 


dynamical 


0.1832(3) 








quenched 


0.2171(1) 


3.42(3) 


0.94(2) 


0.72(1) 



Table 1: Compilation of different parameters (see text) for the dynamical and 
quenched theories at (3 = and [3 = 0.8 on an 8 3 x 16 lattice. 



The corresponding behavior of at /3 = 0.8 for quenched and dynamical 
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fermions is plotted in Figure From the inset of Figure |] it can be seen, that 
as long as the quenched theory is considered the situation is fully compatible with 
eq.(|9|). Thus eq. ([X0|) holds as in the case of (3 = 0. Concerning dynamical fermions, 
again the situation at (3 = 0.8 appears to be in sharp contrast to the (3 = and 
to the quenched could be expected from the properties of the pion norm. 

By approaching the 'critical' value k c {(3) from below with dynamical fermions, 
the ^-dependence of is not linear anymore, i.e. is not compatible with eq.([5|). 
Moreover, close to k c {(3) has a comparatively large finite minimal value which 
would imply that a zero-mass pseudoscalar particle is not contained in the spectrum 
of the theory at this particular coupling. Increasing n beyond k c the pseudoscalar 
mass starts to rise again. Note that in the vicinity of k c the dependence of the 
pseudoscalar mass on k is different for k > k c and k < n c , in accordance with 
the discussion of the pion norm before. 

3 Conclusions and discussion 

We have studied the approach to k c (/3) for two /3-values within the confinement 
phase of the compact lattice QED with Wilson fermions comparing the full the- 
ory with its quenched approximation. We have shown the importance of vacuum 
polarization effects due to dynamical fermions in the context of the chiral limit. 

In the strong coupling limit (3 = the main effect of dynamical fermions seems 
to be a renormalization of the 'critical' value k c , K c yn ^ Kf en . The functional 
dependence of the studied observables on k (respectively m q ) in the limit k — > k c 
compared to the quenched approximation does not change. Our data suggest, that 
at (3 = the pseudoscalar particle becomes massless when k — > k c . 

At (3 = 0.8 the presence of the dynamical ('sea') fermions drastically change 
the transition. There we have found a transition which cannot be associated with 
the zero-mass limit of a pseudoscalar particle anymore, in sharp contrast to the 
quenched case. 

Naively, one would expect that the chiral limit could be established everywhere 
in the confinement phase when approaching the 'critical' line k c ((3). This is not 
the case. Therefore the question arises whether even at (3 = the vanishing pseu- 
doscalar mass can be interpreted as a chiral limit. An alternative scenario for the 
vanishing pseudoscalar mass could be the breakdown of some other symmetry than 
the chiral one. 

Another possible conclusion from these observations is that the calculations - 
analytical and numerical - in the strong coupling limit (3 = can hardly serve for 
the interpretation of the mechanism of the chiral transition at larger /3's (at least, 
in this model). 

These statements need further confirmation, especially on larger lattices. Here 
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work is in progress [10]. It is interesting to what extent the given conclusions can 
be generalized to the case of QCD. 
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Figure captions 

Figure [I]: The pion norm (n) vs. m q at f3 — 0. The curved lines correspond 
to eq.(|) with Co, C\ given in Table [|. 

Figure 0: The unnormalized distribution of n at /3 — 0.8 in the vicinity of k c 
(a) and (p m0 n) in dependence of k at the same value of (3 (b) both for the theory 
with dynamical fermions. 

Figure ^: Counterpart to Figure [I] at (3 = 0.8. 

Figure The K-dependence of m\ for the quenched and dynamical theories at 
(3 = on an 8 3 x 16 lattice. The broken lines represent linear fits. 

Figure |5£ The behavior of m\ around k c at (3 = 0.8. The /t-scale of the smaller 
plot is condensed in order to permit a direct comparison between the dynamical 
and quenched data. The straight line corresponds to a linear fit of the quenched 
data. 
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Figure 1: The pion norm (11} vs. m q at (3 = 0. The curved lines correspond to 
eq.(j8|) with C , C\ given in Table [1]. 
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Figure 2: The unnormalized distribution of II at (3 = 0.8 in the vicinity of k c 
(a) and (p mon ) in dependence of k at the same value of (5 (b) both for the theory 
with dynamical fermions. 
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Figure 3: Counterpart to Figure [T| at (3 — 0.8. 
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Figure 4: The ^-dependence of for the quenched and dynamical theories at 
(3 = on an 8 3 x 16 lattice. The broken lines represent linear fits. 
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Figure 5: The behavior of around k c at /3 — 0.8. The /t-scale of the smaller 
plot is condensed in order to permit a direct comparison between the dynamical 
and quenched data. The straight line corresponds to a linear fit of the quenched 
data. 
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